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We consider a class of percolation models, called Achlioptas processes, discussed in [Science 323, 
1453 (2009)] and [Science 333, 322 (2011)]. For these the evolution of the order parameter (the 
rescaled size of the largest connected component) has been the main focus of research in recent years. 
We show that, in striking contrast to 'classical' models, self-averaging is not a universal feature of 
these new percolation models: there are natural Achlioptas processes whose order parameter has 
random fluctuations that do not disappear in the thermodynamic limit. 
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INTRODUCTION 



THE MODEL 



Percolation is a fundamental problem in statistical 
physics, and the emergence of long range connectivity 
in various percolation models is one of the quintessential 
examples of a phase transition. While many 'classical' 
models and their broad universality classes are nowadays 
well understood [iHl] , in recent years a new class of per- 
colation models has been widely studied. These so-called 
Achlioptas processes 043 are defined via a slight modifi- 
cation of well-studied Erdos-Renyi random graphs, and 
have been of great interest to many physicists [J, IsMSOf 
due to their intriguingly different features. 

One of the most studied properties of Achlioptas 
processes is the evolution of the order parameter (the 
rescaled size of the largest connected component), yield- 
ing several surprises [1, 0, [3, H]- Indeed, certain 
Achlioptas processes were first claimed to have discon- 
tinuous (first order) phase transitions Q , in striking con- 
trast to the typical second order transition observed in 
classical percolation models. Called explosive percolation, 
this was subsequently supported by many researchers 0- 
15l . |21| . but recently it was mathematically rigorously 
proven that the transition is in fact continuous for all 
mean- field Achliopt as p rocesses 0, [l^j , although it can 
be extremely steep |16| . Furthermore, compared to the 
classical cases many of these models seem to behave dif- 
ferently in essential ways, representing a new universality 



class 3 23b J2c ,2£|whose basic features still require fur- 
ther investigatiori[3ll. [s^ . 

In this letter we show that these new percolation mod- 
els can show another surprising behaviour: there arc 
natural Achlioptas processes whose order parameter has 
large random fluctuations, i.e., is not self-averaging. This 
is in contrast to classical models, where the order param- 
eter converges to a nonrandom function in the thermody- 
namic limit. As our nonconvergent examples are from dif- 
ferent universality classes, including the 'explosive' one, 
these also serve as a cautionary tale: when studying a 
wide range of Achlioptas processes, one should not take 
convergence for granted. 



Many 'competitive' percolation models on the com- 
plete graph have been studied, see e.g. 
and the references therein. These Achlioptas processes 
start, as in the classical Erdos-Renyi (ER) model, with 
an empty graph consisting of a large number n of isolated 
vertices, and then sequentially add random edges to the 
graph. In the simplest case in every round two random 
edges ei = viV2 and 62 = V3V4 are picked (rather than 
one), and, using some rule, one of them is chosen and 
added to the evolving graph. Of course, there is nothing 
special about two edges; indeed, Achlioptas processes are 
a subset of the general class of i-vertex rules introduced 
in 0. Here in every round i uniformly random vertices 
vi, . . . ,V£ are chosen, and two of them are connected with 
an edge using some rule (see Q for more variations). 
Note that the resulting percolation processes include the 
ER model, which we obtain by always adding viV2- 

In the following we introduce two different 3-vertex 
rules belonging to different universality classes, always 
writing Ci for the size of the component containing Vi. 
The first rule, the NG rule, proceeds as follows. If all 
three component sizes Ci are equal, add viV2- If exactly 
two component sizes Ci are equal, connect the correspond- 
ing vertices with an edge. Otherwise (if all Ci are differ- 
ent) join the vertices in the two smallest components. 
This rule is named after Nagler and Gutch, who sug- 
gested a slight variant in a different context (personal 
communication). It is a modification of the 'explosive' 
triangle rule introduced in 11, which has a steep 
(but continuous) transition. The rapid growth of the or- 
der parameter in Figure |3] shows that NG also belongs to 
the class of 'explosive' rules. 

The second rule we consider can be viewed as a mod- 
ified ER model and thus we use the shorthand mER. 
Writing Li and L2 for the sizes of the two largest com- 
ponents of the evolving graph, it is defined as follows. 
If the two largest components in the current graph have 
the same size {Li = L2), add viV2- When Li > L2, if at 
least two Ci are equal to Li, connect two corresponding 
vertices; otherwise connect two vertices in components of 



size smaller than Li. Figure [3] indicates that mER be- 
longs to the class of 'nonexplosive' rules, where the order 
parameter evolves rather slowly around the percolation 
threshold (as in the classical ER case). 

RESULTS 

For Achlioptas processes the natural order parameter 

is 

n 

where Li(m) denotes the size of the largest connected 
component after m steps (suppressing in the notation 
the dependence on the rule TZ and on the number of ver- 
tices n) . For classical percolation models it is well known 
that the order parameter converges to a nonrandom func- 
tion in the thermodynamic limit, i.e., that there exists a 
scaling limit p : [0, oo) [0, 1] such that 

\im pn{t) = p{t). (1) 

n— f oo 

Since Li{tn) is random, this actually means that Pnit) = 
Li{tn)/n converges in probability to p{t), see e.g. [7| for a 
more detailed discussion. So ^ asserts: (A) that Pn{t) 
is self-averaging, i.e., that the random fluctuations are 
negligible as n — > oo, and (B) that the expected value of 
Pn it) converges to some deterministic function, i.e., that 
¥,pn{t) does not depend on n in the thermodynamic limit. 

In the following we show that for the NG and mER 
rules the corresponding scaling limits do not exist, i.e., 
that ((H fails. To establish this we estimate the mean 
Hn and standard deviation cr„ of jOn(l) for a range of n 
(based on 10^ runs for each n). Formally we say that 
Pn(t) is self- averaging (at time t = 1) if On/pn — > as 
n — > oo, which is well known to hold for the ER process 
(see also Figure [J). For the NG rule, and especially the 
mER rule. Figure [T] demonstrates that (Jn/ Pn does not 
tend to in the thermodynamic limit; hence these are 
not self-averaging. This is further illustrated in Figure [H 
which also indicates that the expected values of the corre- 
sponding order parameters have a non-trivial dependence 
on n: they seem to oscillate. While for the mER rule the 
amplitudes seem to vanish in the thermodynamic limit, 
for the NG rule we leave it open whether EpJJ^(t) con- 
verges as n — oo or not: if it does, then Figure[2]suggests 
that this only happens for extremely large values of n. 

To put our findings into perspective, note that it is easy 
to construct artificial Achlioptas processes without scal- 
ing limits (using, for example, different rules depending 
on the vertices sampled in the first round). The point is 
that the mER, and in particular the NG rule, are natural 
examples which nevertheless are nonconvergent. Further- 
more, they belong to two very different classes of rules 
(explosive and nonexplosive). This shows that, contrary 
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FIG. 1. (Color onlme). Estimated value of cr„//i„ for pniX) as 
a function of n (each based on 10'' runs) using the NG, mER, 
ER rules and product rule (PR) [J. Self-averaging (ER, PR) 
does not correlate with being 'explosive' (NG, PR). 
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FIG. 2. (Color online). Simulation of the order parameter 
pniX) ~ L\{n)/n as a function of n for the NG and mER 
rules: average and standard deviation (10* runs for each n). 



to classical percolation models, self-averaging is not a 
universal feature of Achlioptas processes. 
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FIG. 3. (Color online). Simulation of the order parameter 
Pn{t) ~ Li{tn)/n as a function of t for the NG and mER 
rules with n = 10^^: five sample runs together with min, max 
and average of 10^ runs. 



ANALYSIS 

We start by collecting some structural properties of 
Achlioptas processes. Remark 9 in Q implies that whp 
(meaning with high probability, i.e., with probability 
tending to 1 as rt — oo) the following two properties hold 
throughout the entire evolution of any ^-vertex Achliop- 
tas process: (i) the vertex set of the union of the £ — 1 
largest components changes by at most o(n) vertices in 
any o{n) steps, and (ii) there exists a function s = s(n) 
such that the ^-th largest component always has size at 
most s = o(n). Although (i) might seem rather tech- 
nical at first sight, it has an important implication: it 
shows that (whp) the size of the largest component can 
'jump', i.e., change by Q{n) vertices in o{n) steps, only if 
there is a step where two linear size components merge to 
form the new largest component (for ^-vertex rules this 
strengthens some of the main conclusions in [31 ) ■ 

Turning to the evolution of the NG and mER rules, 
the key property of both is that they prevent the largest 
component from growing as long as the two largest com- 
ponents have different sizes (Li > L2). But, whenever 
we have two linear size 'giant' components of the same 



size, they merge with constant probability in each step 
(recall that by property (ii) we can have at most two 
such components as here £ = 3). If this happens then Li 
doubles and by (ii) we are left with one giant; all other 
components, in particular L2 , have size at most s = o(n) . 
Now Li is prevented from growing and L2 slowly starts 
growing. After some time L2 emerges to linear size, while 
all smaller components have size at most s = o{n) by (ii). 
When Li and L2 have similar size, they might overtake 
each other repeatedly, but since we have many small com- 
ponents they will merge in o(n) steps, and again Li es- 
sentially doubles by property (i). To summarize, after the 
first linear size giant appears, the remaining evolution of 
Li is essentially governed by discrete doublings. 

The behaviour at the percolation threshold, where the 
first linear size component is formed, seems rather com- 
plicated for both rules. However, the simulation results 
for the NG and mER rules depicted in Figure [3] show that 
the size of the linear size giant component is not deter- 
ministic once it first emerges: note the large spread in the 
order parameter Pn(t) = Li{tn)/n. So, after the perco- 
lation threshold the largest component has random size. 
Since the remaining discrete doublings depend on this ini- 
tial value, this shows that the corresponding Pn^{t) and 
Pn^^{t) exhibit large random fluctuations in their subse- 
quent evolution. Figures [T] and [2] indicate that these do 
not disappear in the thermodynamic limit: the order pa- 
rameters of the NG and mER rules are not self-averaging. 
This 'freezing in' of early variations (the microscopic fluc- 
tuations in the size of the largest component are magni- 
fied and propagated to later in the process) is similar to 
the nonconvergent behaviour of the maximum degree in 
the Barabasi- Albert network model [11], see e.g. |3^ . 

For the mER rule we can formally describe this prop- 
agation to later stages, since its evolution has a close 
connection to Erdos-Renyi random graphs. To see this, 
we first argue that throughout its evolution, up to o{n) 
differences, mER yields a graph of the following form: its 
largest component has size an, where a > 0, and the re- 
maining smaller component sizes follow the distribution 
of an ER graph G(n, X/n) where n = (1 — a)n and A > 0; 
here G{n,p) denotes the binomial random graph model 
in which, starting with an empty graph on n vertices, 
each of the (2) possible edges is included independently 
with probability p (for p — 2m /n^ it is well known to 
be very similar to the uniform ER graph obtained after 
inserting m random edges). The key observation is that 
whenever Li > L2, the mER rule either adds an edge 
with both ends inside the largest component, or one with 
both ends outside. In each case the endpoints are chosen 
uniformly at random. So, between each doubling of Li 
the graph outside the largest component evolves like a 
(rescaled) ER graph. Once Li and L2 have similar linear 
size they merge in o(n) steps, so, for equality up to o{n) 
differences, it sufiices to show that the step in which they 
merge retains the claimed form of the graph. But this fol- 
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lows from the discrete duality principle for G(n, c/n) first 
used in [ssj : the removal of the largest component leaves 
an ER graph G{n', d/n') on n' vertices with d < 1. Sum- 
marizing, the evolution of the mER graph after m = tn 
steps is described by the evolution of the parameters a{t) 
and X{t), which control the size of the largest component 
and the distribution of the smaller component sizes, re- 
spectively. For i < 1/2 we have a = and A = 2t: the 
graph is still close to the ER graph (as all components 
are small). Close to t = 1/2 the parameters evolve in a 
non-deterministic way. Then, as a function of this 'ini- 
tial randomness' at time to > 1/2, the later evolution 
is deterministic and can be described explicitly for any 
t > to- For the mER rule this explains how the lack of 
self- averaging depicted in Figure |3] arises. 

Finally, the 'jumps' of the order parameter in Fig- 
ure [3] show that the NG and mER rules are exam- 
ples of Achlioptas processes which violate (★) continu- 
ity throughout the process, although these satisfy the 
widely studied (f) continuity at the percolation thresh- 
old by 0, [i^l- In fact, since in both rules the largest 
component grows essentially only by discrete doublings, 
it follows that after the percolation threshold both give 
rise to an 'infinite' number of discontinuities in sense (*). 
Property (i) implies that in £- vertex rules such disconti- 
nuities can only arise if two linear size components merge. 
It would be interesting to know if, when restricting to £- 
vertex rules whose decisions (which edge to add) depend 
only on the component sizes (ci, . . . , q), this is the only 
mechanism leading to nonconvergent behaviour. 



CONCLUSIONS 



In summary, we have shown that self-averaging of the 
order parameter is not a universal feature of Achlioptas 
processes: the corresponding scaling limits do not always 
exist. While for any particular rule this can be tested 
numerically via simulations (as in 23, 2^ for example), 
the large interest in such percolation models [1, Isl-fsojl 
indicates that a theoretical investigation of convergence 
in Achlioptas processes is needed, see also [3l|, [s^]- Is 
there a general criterion which allows us to determine 
whether certain rules are convergent or not? Recently 
self-averaging has been rigorously established for a re- 
stricted class of Achlioptas rules [s^: this includes the 
'dCDGM' rule considered in llGl . but not the widely stud- 
ied 0, [i-[ll M, m, EMilillSl product mle (PR). 
Simulations indicate that the scaling limit of the PR rule 
most likely exists (see e.g. Figure [B, but given the sur- 
prises that it has shown so far 0, mi, we can only 
be sure once we have rigorous results. 
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